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SUMMARY 


Equations for analyzing the potential gust-induced overspeed tendency of 
helicopter rotors are presented and discussed. A parametric analysis has also 
been carried out to illustrate the sensitivity of rotor angular acceleration to 
changes in rotor lift, propulsive force, tip speed, and forward velocity. 

The results are in good agreement with a more complex and comprehensive 
nonlinear helicopter simulation analysis and indicate that an increase in 
helicopter gross weight increases the overspeed potential, whereas an increase 
in parasite drag decreases the overspeed potential. 


INTRODUCTION 

The gust response of helicopter rotors has been the subject of study for 
many years and, generally, a rotorcraft has been found to be less sensitive to 
gust- induced airloads than a fixed-wing aircraft (see, e.g., ref. 1>. A recent 
incident, however, raised a question as to the sensitivity to gusts of the 
rotor rotational speed. For the incident considered, the rotorcraft reportedly 
encountered an intense gust which resulted in a rotor overspeed condition. Sub- 
sequent events led to other complications and, ultimately, a forced landing. 

The purpose of the analysis developed and presented herein is to characterize 
the initial gust-induced rotor overspeed tendency. 

In order to examine the potential for a rotor overspeed and to assess the 
sensitivity of the problem, an analysis based on simple blade-element theory 
has been developed. The results of this parametric analysis and the equations 
used are presented in this paper. Correlation of the analysis with a more com- 
plex and comprehensive rotorcraft flight-simulation computer analysis is also 
presented. 


The positive direction of forces, moments, and angles are depicted in 
figure 1 . 

a rotor-blade-section mean lift-curve slope, 5.73/radian 

a] longitudinal flapping coefficient of rotor tip-path plane, radians 

C D ^ R rotor drag coefficient, D R /PTTR 2 (^R ) 2 

C L,R rotor lift coefficient, Lr/PTTR^ (^R) 2 

c P,i 


rotor induced-power coefficient 
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C FfQ rotor profile-power coefficient 

c P,p rotor parasite-power coefficient 

CQ f e engine torque coefficient, Q/pirR^ (£3R) 2 

c q,r rotor torque coefficient, Q/p7TR 3 (ftR) 2 

Cq rotor angular acceleration coefficient, filp/pTTR 3 (fiR) 2 

Dp rotor drag, newtons 

f helicopter equivalent flat-plate drag area, meter 2 

Ip rotor polar moment of inertia, kilograms-me ter 2 

Lr rotor lift, newtons 

Q torque, newton-meters 

R rotor-blade radius, meters 

v rotor mean induced velocity, meters/second 

V airspeed, meters/ second or knots 

Vg gust velocity, meters/second, positive upward 

a angle of attack of rotor control axis, radians 

a Tpp angle of attack of rotor tip-path plane, radians 

<$ rotor-blade-section mean drag coefficient, 0*01 

0 rotor-blade collective pitch at 0.75R, radians 

X inflow ratio, (V sin a - v)/fiR 

M tip-speed ratio, *V/fiR 

p mass density of air, kilograms/meter 3 

a rotor solidity? ratio of total blade area to disk area 

rotor angular velocity, radians/second 
ft rotor angular acceleration, radians/second 2 
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METHOD OF ANALYSIS 


The analysis carried out herein is based on a development and solution of 
the torque balance equations for the rotor and uses the sign convention shown 
in figure 1. The equations are derived and presented in the appendix. The 
basic rotor acceleration equation given in terms of the defined acceleration 
coefficient is 


l 9 c q,r v g 

c « = - C Q,E + C Q,R + ~T— 


( 1 ) 


It is apparent that the rotor can accelerate or overspeed only when the 
net value within the parentheses of equation (1) is negative (i.e., when the 
rotorcraft enters autorotative flight, or when the engine shaft torque supplied 
exceeds the rotor torque required, or when the gust perturbation is large and 
negative in sign). The solution to equation (1) involves solving for the ini- 
tial rotor torque required Cq,r for the appropriate flight state, computa- 

9c q,r 

tion of the torque derivative , and developing and solving an expression 

3a 

for changes in the governor-controlled engine shaft torque supplied c Q r E* 

For the parametric analysis conducted to obtain the trends predicted herein, 
the assumption was made that the shaft torque supplied remained constant at a 
value which equaled the initial torque required; that is, the fuel control 
governor system characteristics allowed the engine to accelerate the rotor when 
a gust disturbance reduced the rotor torque required. This approach has the 
advantages of being simpler to implement and is representative of the actual 
initial response of most governing systems. By using this approach, the rotor 
acceleration potential can be simply expressed as a function of the torque 
derivative and the disturbance gust as 


Cfi - - 


9c Q,R ^g 

3a v 


(2) 


where 


3C 


Q,R 


3a 



C D,R \ 9c L,R 
C L,R J 9a 


^ C L,R 



(3) 


The solution of equation (2) , as developed in the appendix, evolves from 
the expression for rotor torque based on the sum of profile power, induced power, 
and parasite power. This summation is 
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(4) 


a6 <c l ,r) 2 

C Q,R + ^ J + CD ' RU 


(See appendix, eq. (A8) .) With the assumption made herein that the blade- 
section mean drag coefficient 5 is constant, then the torque change due to a 
gust is basically dependent on changes in and C Df R, which are associated 

with the rotor induced and parasite power terms of equation (4) . As illustrated 
in figure 2, these changes are related through the angle a TPP and ^ ts changes. 
The force diagram of figure 2(a) represents a trimmed flight condition. The 
force diagram of figure 2(b) represents changes which may occur when the trim 
is disturbed by the gust velocity Vg. The resultant torque change is thus 
dependent on the relative magnitude and sign of the changes in rotor lift and 
drag. The results of parametric variations in initial trim lift and drag for 
ranges of rotor speed and airspeed are presented in the following section. The 
calculations are not independent of rotor solidity even though the resultant 
acceleration coefficient is normalized by solidity. 


RESULTS AND DISCUSSION 

The analysis developed and presented in the appendix has been used to 
predict potential rotor overspeed tendencies for a range of initial rotor trim 
conditions. As discussed previously, the calculations were based on a constant 
value for engine shaft torque supplied. A sharp-edge, upward vertical gust of 
9.144 m/sec was used as the disturbance velocity V g . Thus, the magnitude of 
the rotor overspeed potential predicted by the analysis is determined by two 
factors: First, the rotor torque sensitivity to changes in angle of attack 

9c q,r a 

; and second, the change in rotor angle of attack due to the gust Aot « V g /V. 

9a 9C Q,R Vg 

As previously indicated, the product — — must be negative for a rotor 

3a v 

overspeed potential to exist. Results from this analysis are presented in 
figure 3 for a realistic range of the parameters C L ,r/ g an< ^ f/^R^ at 
several values of rotor tip speed ^R and airspeed V. For the calculations 
herein, the following characteristics were also assumed: 

a = 0.055 

a = 5.73 

6 = 0.01 


Overspeed Trends 

Figure 3 has been assembled so that the parametric trends are readily 
assessed. For example, trends of the over speed curves with increasing C^r/ 0- 
at constant f/TTR^ are observed from top to bottom of each column of the 
figure, whereas trends of the overspeed with increasing f/^R^ at constant 
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C L r /° are seen in the rows from left to right. Also, within the individual 
plots the effect of changes in airspeed and/or rotor rotational speed are shown 
These trends are reviewed in subsequent paragraphs. 


Effect of Gross Weight 

A comparison of the columned plots of figure 3 (e.g., figs. 3(a), 3(d), 

3(g) , and 3 ( j ) ) shows that, for a given helicopter with fixed flat-plate drag f, 
an increase in gross weight, as reflected by C L f R/ a , increases the maximum 
rotor overspeed potential due to the gust. The airspeed at which the maximum 
over speed potential occurs also increases with C L ^ R /a. An assessment of equa- 
tion (3) , in consideration of these trends, indicates that the dominant factor 
C L,R* The first term C L ^ R /y arises from the induced-power term of equa- 
tion (4) and, consequently, increases in gross weight are reflected in an 
increase in the induced-power contribution to the torque derivative. For an 
overspeed to occur under the conditions cited, the torque derivative must be 
negative in sign and the induced term must not have a negative value. Thus, 
for conditions shown in figure 3, where an over speed potential is predicted, 

the predominant term is the third term in equation (3), that is, yCL fR ll + 
which arises from the parasite-power term of equation (4) . \ 



Effect of Flat-Plate Drag 

The row plots of figure 3 (e.g., figs. 3(a), 3(b), and 3(c)) illustrate 
that for a helicopter with constant gross weight, an increase in flat-plate drag 
area f/TTR^ decreases the maximum over speed potential at any given airspeed. 
This decrease results from the second term of equation (3) which, by combining 

C D,R f y 2 

with = - — (see appendix, eq. (A10)), can be shown to be propor- 

a ttr 2 2a 

tional to fy^ and positive in sign. Thus, increases in the trim flat-plate 
drag decrease the overspeed tendency. Equally important to the overspeed trend 
is the tip-speed-ratio term y^, which causes large changes with speed. In par- 
ticular, at low-lift and high-drag conditions, as shown in figure 3(c), the 
increase in magnitude of the second term of equation (3) with speed is predom- 
inant in controlling the overspeed trend. However, at the high-lift and low- 
drag conditions illustrated in figure 3(j), the importance of the second term 
is diminished and the overspeed trends are again dominated by the third term 

/ a ai \ 

yCT r 1 + I of equation (3). In any case, however, the airspeed at which 

' \ / 

the maximum overspeed potential occurs is inversely affected by the flat-plate 
drag . 


Airspeed for Maximum Potential Overspeed 

As previously mentioned, the airspeed at which the maximum rotor over speed 
potential exists varies directly with the helicopter gross weight and inversely 
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with the flat-plate drag. For the range of conditions covered in figure 3, the 
curves show that the maximum overspeed potential occurs anywhere from 55 knots 
at the lowest gross weight and flat-plate drag condition of figure 3(a) to 
slightly over 140 knots at the highest gross weight and lowest drag condition 
of figure 3(j). This wide range in speed, however, is dictated by the rela- 
tively large range of lift-drag conditions evaluated, and the more represen- 
tative conditions are those typified by the midrange values (i.e., figs. 3(e) 
and 3(h)). For these conditions, the airspeed for maximum rotor overspeed 
potential occurs between 60 and 100 knots. This speed range is of particular 
significance because a number of military helicopter operation manuals specify 
a gusty or turbulent air penetration velocity of 80 knots. This limit is dic- 
tated by other considerations, but for those helicopters which are susceptible 
to rotor over speed and, in particular, those with low rotor polar moment of 
inertia and sluggish governor response characteristics, the rotor-speed fluc- 
tuation may be disconcerting. 

To illustrate the significance of the overspeed parameter, consider the 
configuration 


c L,R/ a = °* 08 

f/ ttR 2 = 0.009 

V = 80 knots 

ftR = 198.12 m/sec 

R = 4.01 m 

a = 0.055 

Ip = 338.86 kg-m 2 


From figure 3(h), C^/a = 2.25 x 10 -8 and therefore 


. pirR 2 (ftR) 2 RCF . 

Q 3.56 rad/sec^ 

CT Ip 


Since the initial rotational speed is 49.39 rad/sec , the predicted over speed 
tendency represents a potential 7.2-percent increase in rotor angular velocity 
for each second the acceleration persists. 
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Comparison of Simple Theory With Nonlinear Helicopter Theory 

In order to validate the trends established by using the simple analysis 
developed herein, a comparison has been made between these simple blade-element- 
theory results and results obtained from a nonlinear helicopter simulation 
analysis described in reference 2. The nonlinear results were obtained by 
using the maneuver option of the program and allowing the simulated vehicle to 
encounter a 9.144-m/sec, sharp-edge, upward vertical gust. In order to simu- 
late sluggish governor response characteristics, the value of the engine torque 
supplied to the main rotor was fixed throughout the maneuver at the value cal- 
culated as that torque required for the initial trimmed level flight. The 
rotor angular acceleration was determined from changes in the rotor rotational 
speed between the time the rotor first encountered the gust and the time at 
which the rotor had penetrated the gust a distance of two rotor diameters. This 
penetration interval provided a common basis for the calculation of Cft by 
establishing a point in the maneuver where the rate of change of rotor rota- 
tional speed had essentially stabilized. 

Two comparisons are made. Figure 4 shows a comparison of the simple 
results with the nonlinear predictions using steady, two-dimensional blade- 
section aerodynamic characteristics in the rotor analysis. The over speed trends 
predicted by both methods are generally similar, although the nonlinear analysis 
indicates a more dramatic reduction in overspeed potential at the higher air- 
speeds. Because the nonlinear analysis accounts for the effects of rotor blade 
stall and compressibility and the present analysis does not, the torque deriva- 
tive of the nonlinear analysis is probably more susceptible to increases in 
torque due to stall and thus does not show the same tendency to overspeed at 
high forward speeds. 

As previously mentioned, the nonlinear analysis uses steady aerodynamics 
in computing rotor performance, and thus the rotor derivatives may be affected 
by the stall characteristics. In order to assess this effect, the nonlinear 
analysis was rerun by using an unsteady aerodynamics option (identified as BUNS 
in ref. 2). This unsteady option, which is probably a more realistic represen- 
tation of transient rotor aerodynamics, results in a delay in blade-section 
lift stall and a retardation of section drag increase. The results of the 
nonlinear analysis with the unsteady option are also shown in figure 4. Inclu- 
sion of the unsteady aerodynamics effects generally results in better overall 
correlation with the linear analysis. This is perhaps not too surprising 
because the unsteady aerodynamic method used in reference 2 tends to maintain 
a more constant lift-curve slope and to produce less drag growth as rotor-blade- 
section angle is increased. Both characteristics are inherent in simple blade- 
element theory, and therefore the present analysis provides a very reasonable 
trending of the over speed tendency of a helicopter rotor. 


CONCLUSIONS 

The simple analysis of the gust-induced overspeed trends of helicopter 
rotors reported herein has identified the basic parameters affecting rotor 
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overspeed and has provided a range of trending curves to illustrate the more 
important parameters* Specifically, analysis of the results leads to the 
following conclusions: 

1 * An increase in helicopter gross weight increases the overspeed 
potential* 

2* An increase in helicopter parasite drag decreases the overspeed 
potential* 

3. The airspeed at which the maximum overspeed potential occurs varies 
directly with gross weight and inversely with flat-plate drag area. 

4* For representative values of helicopter lift and drag characteristics, 
the maximum overspeed potential occurs in the 60- to 100-knot airspeed range, 
which brackets the generally recommended 80-knot gust penetration airspeed 
limit* 


5* The correlation of the simple analysis with a more complex and costly 
nonlinear helicopter simulation analysis, which included an unsteady aerodynam- 
ics routine, indicated that the simple analysis trends are very representative* 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
May 24, 1978 



APPENDIX 


DEVELOPMENT OF EQUATIONS FOR ROTOR-OVERSPEED ANALYSIS 

The basic equation necessary to determine the over speed potential of a 
rotor system is based on the dynamical equation of the rotor shaft torque 
balance £q + I p n = 0 expressed as 


. p7TR 3 (QR) 2 / 

Q = — I Cr 

Ip \ ^ 


3Cr 


(AT) 


where Cq^ r is the initial trimmed value of rotor torque required. An accel- 
eration coefficient is defined as 


flli 


CQ = 


PTTR 3 (ftR) 2 


(A2) 


Equation (AT) may then be written as 


CQ = - 



+ c 


Q'R 


+ 


3CQ,R Vg 
3a v 


(A3) 


The development of the equations necessary for solving equation (A3) are based 
on the simple blade-element theory as set forth in reference 3. 


Rotor Torque Coefficient 

Because the rotor shaft torque coefficient and the shaft power coefficient 
are equivalent, the expression for the rotor torque coefficient may be written 
as 


C Q,R “ c P,o + c P,i + c P,p (ref. 4) 


(A4) 
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where, by integration, equation (26) of reference 5 may be expressed as 


c P ,o - -n + 3 ^ 2 > 


(A5) 


With the assumption that the rotor thrust coefficient and rotor lift coefficient 
C L R are equivalent and that [l + (A/y) 2 ] 1 / 2 « l, equation (7) of reference 4 
may be expressed as 


c P,i 


c l,r' 

2y 


(A6) 


and with the additional assumption that cos a » 1, equation (8) of reference 4, 
may be expressed as 


C 


PfP 


f y 3 

HR 2 2 


-c d ,r^ 


(A7) 


Substitution of equations (A5) , (A6) , 


'Q,R 


a<$ 2 

= — n + 3y 2 ) + 


2y 


and (A7) into equation (A4) yields 
" C D,R^ 


(A8) 


or 


C Q,R 


a 


6 / c l,r\ 2 a C DfR 

-(1 + 3y 2 ) + y 

8 \ a J 2\l o 


(A9) 


By expressing C DfR /tf as 


C D,R f y 2 

a TfR 2 20 


(A10) 
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equation (A9) becomes 


C Q,R 


a 



+ 3y 2 ) 



f y 3 
TrR 2 2a 


(All) 


Rotor Torque Derivative 

The rotor torque derivative is obtained by a partial differentiation of 
equation (A8) with respect to rotor angle of attack. The resulting expression 
is 


9C Q,R c L,R 9 Cl,R 9c d,r 
9a y 9a 9a 


(A12) 


With the assumption that the resultant force vector is perpendicular to 
the rotor tip-path plane (see fig. 2), then 


C D,R 

C D,R * Ci^rOtPP or a TPP " (A13) 

C L,R 


Thus, 


9c d,r 

8a 


9c l,r 

8a 


a T pp + c l,r 


9a TPP 


8a 


(A14) 


Since, from figure 1 


a T pp = a + a ] 


(AT 5) 


then 


9aTPP 

9a 


9a*j 




9a 


(AT 6) 


n 
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APPENDIX 

Equations (A12), (A13) f (A14), and (A16) can be combined to give 


9c q,r / c l,r c d,r\ 9c l,j 


3a 


M y C LfR / 9a 


/ d &1 

' yCl " R 1 + 5 kT 


(A17) 


Rotor Lift Derivative 

The rotor lift derivative may be obtained from 


ca 

C L,R = ~ 


0^- + y 2 ) + X 


(ref. 3) 


(A18) 


and 


X 


V sin a 

V 

V sin a 

0.R 

fiR 

£2r 


C L,R 



(AT 9) 


or 


X = pa - 



(A20) 


By substituting equation (A20) into equation (A18) to solve for C LfR 
and by taking the partial derivative with respect to a, the lift derivative 
may be expressed as 


9c L,R 2\ i 2 


da 8y 

1 + — 
aa 


(A21) 
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The rotor 


Rotor Flapping Derivative 
flapping derivative may be obtained from 


a l 




(ref. 3) 


(A22) 


Combining equations (A20) and (A22) and differentiating and substituting equa- 
tion (A21 ) yields 



(A23) 


Overspeed Solution 

Equations (All), (A17), (A21), and (A23) are used in the solution of equa- 
tion (A3) for the overspeed coefficient. The computation requirements are min- 
imal and can be carried out on a programmable pocket calculator. For the com- 
putations of this paper the engine shaft torque was assumed to be a constant 
value equal in magnitude but of opposite sign to the initial rotor torque 
required. 
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Figure 1 Notation showing positive direction of forces and angles. 



-c 


D,R 

I C 




■«TPP \ 


V 


L,R 


Drag = (l/2pV 2 )f 
Gross weight 


(a) Trimmed force diagram. 



(b) Disturbed force diagram. 

Figure 2.- Graphical representation of trimmed and disturbed 

rotor force components . 
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